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ABSTRACT 

Central peaks in the iron abundance of intracluster plasma are a common feature of cooling- 
core galaxy clusters. Although centrally localized, these abundance peaks have a much broader 
profile than the stars of the central brightest cluster galaxy (BCG), which produce the excess 
iron, indicating that metal-enriched plasma is transported out of the BCG by some process such 
as turbulent diffusion. At the same time, cooling-core clusters are likely heated by central active 
galactic nuclei (AGNs) by means of jets, cosmic-ray bubbles, and/or convection. The recent 
AGN-driven convection model of Chandran & Rasera predicts the turbulent velocity profile in 
a steady-state cluster in which radiative cooling is balanced by heating from a combination of 
AGN-driven convection and thermal conduction. We use the velocity profiles from this model as 
input into an advection/diffusion model for the transport of metals in the intracluster medium, 
taking the iron to be injected by the BCG. We compare the results of our model to XMM and 
Chandra observations of eight clusters. Assuming a constant turbulence level over a cluster's 
lifetime, the turbulent velocities in the model can explain the observed abundance profiles in 
only five of the eight clusters. However, we go on to develop an analytic fit of the turbulent 
velocity profile as a function of the AGN power. We then deduce for each cluster the average 
AGN power (during the past ^10 Gyr) required to match the abundance profiles. The required 
average values are between 10^^ and 2 x 10^^ erg.s"^, while the present AGN powers span a 
much larger range from 6 x 10*"'^ (Virgo) to 2 x 10^^ erg.s""'^ (Hydra A). Our results suggest that 
AGN-driven convection can account for the observed abundance profiles if the AGN power varies 
over a cluster's lifetime between Virgo-like and Hydra-A-like values, with average values in the 
above-quoted range. 

Subject headings: galaxies: abundances — diffusion — convection — galaxies: clusters: general — cooling 
flows — galaxies: elliptical and lenticular, cD — methods: analytical - galaxies: clusters: individual 
(Perseus, Hydra A, Sersic 159-03, Abell 262, AbeU 1795, Virgo, AbeU 496, Abell 4059) 



1. Introduction 

In many clusters of galaxies, the radiative cool- 
ing time at the cluster's cente r is m uch shorter 
than the cluster's age (Fabianl 119941 ). Neverthe- 
less, high-spectral-resolution X-ray observations 
show that ver y little plasm a actually cools t o low 
temp e ratures (iBohringer et al.. .2001; .David et al 



2001 



2001 



Molendi fc Pizzolatol 2001: IPeterson et al 



2003'; 'Ta mura et all l2001bt iBlanton et al 



20031 ). This finding, some times referred to as 



the "cooling-flow problem," strongly suggests that 
plasma heating approximately balances radiative 
cooling in cluster cores. 

Although a number of different heat sources 
have been considered in the literature, there is 
growing interest in the role of central active galac- 
tic nuclei (AGNs). The importance of AGN heat- 
ing or "AGN feedback" is suggested by the ob- 
servation that almost all clusters with strongly 
cooling cores possess act i ve central radio so urces 
(jBurnsI llQQOl : iBall et al.l Il993t lEilekl \2004 ) and 
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by the correlation between the X-ray luminos- 
ity from within a cluster's cooling radius and 
the mechanical luminosity of a cluster's central 
AGN (|Birzan et al l booil: lEilekl 12004 . One of 
the main unsolved problems regarding AGN feed- 
back, however, is to understand how AGN power 
is transferred to the diffuse ambient plasma. A 
number of mechanisms have been investigated, in- 
cluding Compton heati ng ( Binncv & Tabor 1995 



Ciotti fc Ostrikeij ll997L I2OOII: ICiotti et al.l I200I 



Sazonov et al 
1993; 



20051). 



shocks ( Tabor fc BinnevI 
magnetohydrody- 



B innev fc Taboil I1995D . 
nami c wave-mediated plasma heating by cosmic 



rays te oehringer fc Morfilll ll988l : lRosner fc Tuckeij 
[l989; Loewenstein et al. Il991 ). and cosmic-ray 



bubb l es pro duced by the central AGN (IChurazov et al 
2001', '2O O2I: iRevnolds et all I2OO2I : iBruggenI l20oi 



Reynol ds et al. 20051 ) . which can heat intracluster 

plasma by generating turbul e nce ([Loewenstein fc Fabian 



1990 : IChurazov et all 12004 ICattaneo fc Tevssieil 



2007) and sound wav es ( Fabian et al. 2003at 



■Ruszkowski et al.ll2004allbl) and by doing pdV work 
jBegelman' '20011: 'R uszkowski fc Begelmanl 12002 : 
[Hoeft fc Briiggcn 200^^ 

Another way in which central AGNs may 
heat the intracluster medium (ICM) is by ac- 
celerating cosmic rays that mix with the intr- 
acluster plasma and cause the ICM to become 
convectively unstable. A steady-state, spheri- 
cally symmetri c model based on t his idea was 
developed by IChandran fc Maron (|2p04) and 
subsequently refined by 'Chandrad (j2005l) and 
[Chandran fc Rasera (2007). In this model, a cen- 
tral supermassive black hole accretes hot intra- 
cluster plasma at the Bondi rate, and converts a 
small fraction of the accreted rest-mass energy into 
cosmic rays that are accelerated by shocks within 
some distance rgourco of the center of the cluster. 
The resulting cos mic-ray pressure gradient leads 
to convection (see Chandran fc Dennii ( 2006l l and 



Dennis fc ChandranI (120081 )). which in turn heats 



the thermal plasma in the cluster core by advect- 
ing internal energy inwards and allowing the cos- 
mic rays to do pdV work on the thermal plasma. 
The model also includes thermal conduction and 
cosmic-ray diffusion (viscous dissipation turns out 
to be smaller than the other forms of convective 
heating at all radii) and assumes a steady state 
in which the net heating rate balances radiative 
cooling throughout the cluster. The model uses 



mixing-length theory to describe convection and 
its effects on the ICM and predicts a self-consistent 
profile for the rms amplitude of the turbulent ve- 
locity. By adj usting a single parameter in the 
model (rsource), IChandran fc Raseral ( 20071 ) were 
able to achieve a good match to the observed den- 
sity and temperature profiles in a sample of eight 
clusters. In several of the clusters in this sample, 
compact cooling flows arise within the central few 
kpc of the clusters because the rate of radiative 
cooling peaks much more sharply near the clus- 
ter center than either the convective or conduc- 
tive heating rates. At even smaller radii in these 
clusters, the cooling flow makes a transition to a 
Bo ndi flow, in a manne r similar to that described 
by iQuataert fc NaravanI (|2000l ). The size of the 
central cooling flow plays a role in regulating the 
mass accretion rate within the central Bond i flow, 
as described bv IChandran fc Raseral (^2007^ . 

In this paper, we explore the connection be- 
tween this AGN-driven-convection model and the 
observed properties of the iron abundance pro- 
files of cooling-core clusters. While observations 
of non-cooling-flow clusters show a nearly constant 
abundance profile, cooling- core clusters show very 
peaked iron distributions (jPe Grandi fc Molendi 
i2001i) . Observations of the relative abundances of 
oxygen, silicon and iron suggest that the produc- 
tion of iron in these abundance excesses is domi- 
nated by SNIa (and possibly winds) f rom the cen- 
tral brightest cluster galaxy (BCG) (lEttori et"al 



2OO2I: iMatsushita et alJl2002Hde Plaa et al.ll2006l ). 



If cooling cores are preserved over a tim escale 



longer than 5 Gyr (|B5hringer et al.ll2004l ). then 



the observed amounts of excess iron (of order 
10^ Mq) within the cluster core are compatible 
with t he amounts produced w ithin the BCG by 
SNIa (ICappellaro et all [l999l) and stellar winds 



( Ciotti et al.l 1991 ). However, the shape of these 



abundance profiles is still a mystery. The distri- 
bution of iron is much broader than the distribu- 
tion of the stars that produce the metals, which 
indicates some additional processes are needed to 
transport the metals out of the BCG into the sur- 
rounding ICM (,Rebusco et al.ll2005l ). 



Recently, iRebusco et al.l (|2005l ) developed an 
analytical model of metal injection by SNIa and 
diffusion by turbulent gas motions. They sug- 
gested that the dissipation of the same stochastic 
gas motions would produce the heating required 
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to sol ve the coolinK-flow problem Using this 
model, Rebusco et al. ( 20051 2006h : Graham et al 



( 20061 ) found that diffusion coefficients of order 
X028 _ 2^q29 cm^.s"^ are required to match the 
abundance profiles of cooling-core clusters, while 
spatial scales ~ 10 kpc and velocities of the order 
of few 100 km.s~^ are needed to compensate the 
gas cooling. 

Although these studies offered an explanation 
for the observed abundance profiles as well as a 
solution to the cooling-flow problem, they did not 
explain how the rms amplitude of the turbulent 
velocity, u, is determined, or how the dependence 
of u on the radial coordinate r is determined. On 
the other hand, the AGN-driven-convection model 



of IChandran fc Raseral (|2007[ ) provides a physics- 
based theoretical framework for calculating u{r). 
In this paper, we combine this model for intra- 
cluster turbulence with the m odel of metal injec- 
tion of i Rebusco et all (|2005[ ) in order to under- 
stand the production and transport of metals in 
the ICM. We describe these models further in sec- 
tions [2] and [3l In sections [4] and [5] we apply these 
models to a sample of eight clusters. We discuss 
the possible implications of our results for the vari- 
ability of the AGN power during a cluster's life- 
time in section |6l and summarize our conclusions 
in section [71 

2. Model of iron injection by SNIa and 
winds 

From the relative abundance of O, Si and 



Fe, Ettori et a 



( 20021 ): [Matsushita et al.l (|2002D : 



de Plaa et aLr(|2006h have shown that SNIa dom- 



inate the iron enrichment in these abundance 
peaks. After removing the contribution from SNII 
thought to originate from the early formation of 
the BCG as well as the background abundance 
from other galaxies, the remaining central iron 
excess originates mainly from SNIa and stellar 
wi nds of the central BC G . Following th e work 
of iBohringer et al.l (|2004l) : iRebusco et all (|2005l 
20061 ). we estimate the contribution of SNIa as. 



dp 



Fe 



dt 



= 10 



-12 f sr 



SNIa 



SNU 

fc 



VFe \ / PL 



— ) MQ.yr- 



.kpc 



-3 



(1) 



with sr the present supernova rate in SNU (super- 
novae per century and per 10"^° Lq), rjpe — 0.7 Mq 
the iron yield per SNIa, pL the blue luminosity 
density, t the cosmic time, tn — 13.45 Gyr the 
current age of the universe and k the exponent 
which describes how the supernovae rate increased 
in the past. The w i nd contribution is tak en from 
Ciotti et all (|l99ll ): IRebusco et all (|2005h . 
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Mo.yr-i.L; 




PL 



L|.kpc-3 



— j Mo.yr-i.kpc- 



with = 2.8 X 10~3 the mean iron mass frac- 
tion in the stellar winds of an evolved population, 
"rn-uiind = 2.5 x 10^-^^ MQ.yr^-^.Lg^"^ the present 
star mass loss rate of a ~ 10 Gyr old population 
and, a = 1.3 specifying the time evolution of the 
star mass loss rate. The hidden parameter here is 



tage, which is the age of the BCG. Eq[T]and EqO 
are indeed only valid from tn ~ tage to tu- Be- 
fore the formation of the galaxy, the production 
of met als is of course as sumed to be null. We as- 
, Hernauis"3 ( 1990f ) profile typical of elliptical 



sume a 



galaxies. 



PL{r) 



27r r[a -t- rY ' 



(3) 



where a — re/1.8153 and is the effective radius 
containing half of the projected luminosity. 

We constrain the parameters sr, k and tage so 
that the total observed amount of iron Mpe equals 
to the total amount of iron from the model inside 
the radius corresponding to the last bin of the 
abundance observations (located at a radius rt). 
Using this normalization procedure makes the re- 
sulting abundance profile quite insensitive to the 
particular values of our poorly known parame- 
ters because the only important ingredient for our 
purpose is the light profile. For example, varying 
k between and 2 changes the final abundance 
by only a few perce nt (with the d if fusion coeffi - 



cients suggested bv [ Rebusco et al. ([20051 



Graham et al 




Renzini et al.l 



2006) 



19931 ): 



Rebusco et al.[ ([2006D suggest a value of k be- 
tween 1 and 2, we pick a value of A: = a (that 
is the time dependence SNIa iron injection and 
wind iron injection are the same). In this way, we 
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factorize the time dependence and adopt a new 
simple expression, 



dt 



10 



-12 



VSNUJ 



VFe 

Mo 



PL 



Lg.kpc- 



Mo.yr ^.kpc 



(4) 



with a = 1.3 and sr^ff an effective supernovae rate 
which includes the wind contribution. This effec- 
tive supernova rate will be determined directly 
from the observed mass of iron and is no longer 
a free parameter. The advantage of using Eq|4]is 
that our problem of metal production and trans- 
port depend now linearly on sr^ff. We therefore 
don't need to run our solver for the transport equa- 
tion multiple times to get the right sr. Instead, 
we run our solver once and then adjust the value 
of srcff afterwards (thereby multiplying the abun- 
dance profile by a constant) to match the total 
observed iron mass. 

The last parameter is the age of the BCG, which 
affects our results to some degree since an older 
galaxy allows more time for iron production and 
diffusion. A typical galaxy age is given by the stel- 
lar population of the BCG of order ~ 8 — 12 Gyr 
( Jimenez et aLllioOTl ). We adopt a fixed age in this 
range for all our clusters of 9 Gyr. The metal en- 
richment model is now entirely determined by the 
observational data Lb, r^., Mp^. 

3. Diffusion of metals by AGN driven con- 
vection 

We assume that there is steady-state convec- 
tion in the ICM as d escribed by the AGN-drive n 
convection model of IChandran fc Raseral (|2007l ). 
Metals in the ICM are then advected as a pas- 
sive scalar by the turbulent fiow. The advection 
equation for the transport of metals is given by, 



dpa 



-V.(pva) 



dt 



(5) 



with p the gas mass density, a the mass fraction of 
iron, and v the velocity field given by the AGN- 
driven convection model. Using mixing length the- 
ory and assuming spherical symmetry and a sta- 
tistical steady state, we average equation ([5]) to 



obtain the equation 
d < p X a > 



+ — — 

J.2 Qj, 



dt inr'^ dr 

2t^ d <a> 

r^D <p> — 

or 



M d < a> dpFe 

- < > 



dt 



(6) 



where (. . .) indicates an ensemble average, M is 
the mass accretion rate which is independent of 
the radius and is equal to the Bondi rate near the 
center (at r « 0.2 kpc), < %^ >« ^ is the 
above metal production rate from EqHl < p >— p 
is the gas density, < a >= a is the mass fraction 
of iron, 

D = 0.5luNL (7) 

is the diffusion coefficient, I = 0.4r is the mixing 
length, and unl is the radial component of the 
turbulent velocity from the mixing-length-theory- 
based A GN-driven convectio i i mod el, as calcu- 
lated by IChandran fc Rasera ( 2007 ). The factors 
0.5 and 0.4 ar e choosen to be cons i stent with the 
value used bv IChandran fc Rasera (2007) for the 
different energy fluxes. The factor of 0.5 in equa- 
tion ([7]) is much larg er than the numerica l coeff i- 
cient in Equation 8 o nnis fc ChandranI (|2005l ). 
which reads I?cddy = 0.11/71/^. However, in this 
last equation, U is the rms value of the full ve- 
locity vector, not just its radial component. Sim- 
ilarly, L is the full correlation length of the tur- 
bulence, whereas the mixing length I is just the 
radial component of the typical displacement of 
a fluid clement before it is mixed into the sur- 
rounding fluid. (The quantity ^ accounts for the 
possible reduction in Deddy if the motions are os- 
cillatory and only weakly turbulent.) If one sets 
u ~ U/VS and I ~ L/V3, then D = 0.167 LU 
in equation (O, which is sim ilar to equation 8 



of Dennis fc ChandranI (120051 ). It should also be 
noted that if we were to choose a smaller nu- 
merical constant in equation (O and throughout 
the m ixing-length model of IChandran fc Rasera 
(|2007l ). the resulting turbulent velocities in the 
model would become larger so that the energy 
fluxes and convective heating were still able to 
compensate for radiative cooling. The resulting 
value of D would not be greatly altered, since it 
is closely tied to the rate of convective heating, 
which is fixed in the AGN-convection model by 
the requirement that heating balance cooling in 
steady state. 
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The evolution of the abundance depends on 
two terms. The first term is related to the ac- 
cretion rate M and corresponds to an inward 
advection of the abundance profile at a speed 
Vj^{r) = M / {Anr^ p) . If we were to assume a large 
accretion rate as in the cooling-flow model (e.g., 
M = 100 - 1000 Mo.yr-i) the accretion of low- 
abundance gas from larger radiu s would destroy 



the central abundance excess (jBohringer et al 
20041 ). However, in the AGN-driven convection 
model, M < lM0yr~^, and the effect of in- 
flow is small except within the central few kpc, 
where the l/{pr^) dependence of v^j{r) causes 
Vj;^ to become large. The second term corre- 
sponds to diffusion resulting from the random 
turbulent velocity field. This term dominates 
and smoothes the abundance profile on a scale 
r « 58 X [(D/lO^^cm^s-i) x (tage/lOGyr)]^-^ kpc, 
which is comparable to the typical half-mass ra- 
dius of the central iron excesses. 

We have solved Eq[6] taking a null initial abun- 
dance and assuming a null second derivative of 
the iron abundance profile for the boundary con- 
ditions. The discretization of the advection part 
doesn't require any artificial diffusion since the 
physical diffusion dominates for our case. We have 
run our code to compute the iron abundance pro- 
file of 8 cooling-core clusters. 

4. Properties of the sample of 8 cooling- 
core clusters 

Our sample consists of 8 coohng-core clus- 
ters for which we have the turbulent velocity 
profile from the AGN-driv en convection model 
(jChandran fc Raseral l2007l ): Virgo, Abell 262, 
Sersic 159-03, Abell 4059, Hydra A, Abell 496, 
Abell 1795 and Perseus. The redshifts of these 



clusters are taken from lKaastra et al.l (j2004[ ). The 
properties of these clusters and their brightest 
cluster galaxies (BCGs) are summarized in table 
[T] and table O We use the current ACDM cos- 
mology (iJo = 70 km.s^^.Mpc"^, ri„i ~ 0.3 and 
VLa = 0.7). 

4.1. Photometric properties 

The photometric properties we are interested in 
are the total luminosity L^, and the effective ra- 
dius fe, which contains half of the projected flux. 
We get the name of the brightest cluster galaxies 



(BCG) in each cluster from lRaffertv et all (|2006l ) 

and the NED database, and then compute the to- 
tal blue luminosity using the Hyperleda database 
(jPaturel et al.l |2003j) proceeding as follows. We 
take the apparent corrected blue magnitude of the 
BCGs, compute the luminosity distance, and the 
associated distance modulus from our redshift and 
the ACDM cosmology. We finally compute the 
blue absolute magnitude, and the blue luminosity 
assuming Mq = 5.47. 

Finding the effective radius i s more challenging 



r mdmg the enective radms i s more cnailengmg . 
We use three different sources: ISchombertl ( 19871) 



(Perse us, Abell 1795, and Virgo). Icraham et al 



(|l996l ) (Abell 262 and AbcU 496) and Hyperleda 
(Hydra A, Sersic 159-03, and Abell 4059) In th ese 
references an R^/^ law ( de Vaucouleurd 1953f ) is 
adjusted to the observed profile from which the 
best fit effective radius (or the effective surface 
brightness in Hyperleda) is deduced. We also con- 
vert to our cosmology using the ratio of angular 
diameter distance for our cosmology and redshift 
to the angular diameter distance for the authors' 
cosmology and redshifts. The results for and 
Lcr are shown in table [T] 

We note that we neglect the co res observed 



on small scales ( Carollo et al. 1997[). the de via- 



1996h 



on 



tion from the YO'/'^ law ( Graham et al 
intermediate sc ale, the large en velopes observed 
on large scales ( Schombert 1988 h and the elliptic- 
ity of the profiles. We do not hide here that the 
uncertainties on and Lb are important. For- 
tunately, the result on the shape doesn't depend 
on Lb (due to our renormalization procedure) and 
has only a little dependency on r^, for most of the 
clusters because the diffusion scale is larger than 



4.2. Gas density and turbulent velocity 
profile 

We take the eight clusters in our sample to 
have the density profiles and (as discussed above) 
turbulent-velocity profiles obtai ned from the 
AGN- driven convection model of .Chandran fc Raseral 
(120071 ). The model densi ty profiles are very s imi- 



lar to the observations of ^Kaastra et al. 



as 



shown in Figure 1 of Chandran fc Raseral ( 2007 ). 
We note that the AGN-driven-convection-model 
solutions were obtained by varying the size of the 
cosmic-ray acceleration region in the model in or- 
der to achieve a best-fit to the observed density 
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Fig. 1. — Diffusion coefficient (due to the stochas- 
tic turbulent motions of the gas) for our 8 clus- 
ters. Thi s conies from the AGN-driv en convection 
model of Chandran fc Rasera 



and temperature profiles of the eight clusters in 
our sample. The model solutions were therefore 
not fine-tuned to solve the independent problem 
we are investigating now; the shapes of the abun- 
dance profiles. From the turbulent radial velocity 
profile unl, we deduce the diffusion coefficient of 
EqEl D(r) = 0.2r x unl- 

The profiles of the diffusion coefficients are pre- 
sented in FiglH The diffusion coefficients peak 
around a radius r^ax where the diffusion reaches 
its maximum Dmax- We divide the clusters into 
4 levels of diffusion. Virgo and Abell 496 have 
the smallest diffusion coefficients, with Dmax ~ 
8 X lO^^cm^ .s~^ . They are the "very weak feedback 
clusters". AbeU 4059 and AbeU 262 have small 
diffusion coefficients {Dmax ~ 5 x lO^^cm^.s"^) 
and are called "weak-feedback clusters" . Perseus 
and Abell 1795 have intermediate diffusion co- 
efficients {Dmax ~ lO^^cm^.s"^) and are called 
"moderate feedback clusters". Finally, Hydra A 
and Sersic 159-03 have large diffusion coefficients 
{Dmax ~ 2.5 X lO^^cm^.s"^) and are the "strong- 
feedback clusters". We summarize the value of 
Dmax and Vmax in table El 

4.3. Observed abundance profiles 

We find the observed abundance profiles in the 
literature (helped in this way by the BAX clus- 



ter database ( Sadat et al. 20041 )). All the data 
points come either from Chandra or XMM ob- 
servations. We converted all the radii to kpc us- 
ing the angular diameter distance from KCDM 
cosmology. Finally, we converted all the dif- 
fere nt solar abundance s to th e value published 
by lAnders k. Grevess3 ( 1989( ) (where the so- 



lar abundance relative to H is 4.68 x 10 ^ in 
number). The references used are the foUow- 
ing: Virgo (iMatsushita et all |2002|). Abell 262 
(IVikhlinin et al.ll2005l). Sersic 159-03 (Ide Plaa et al. 
;2006h. Abeir4059 (laioi_etalJ l2004h. Hvdra A 
(iDavid et aP l200lh . Abell 496 (iTamura etU 



2001ah. Abell 1795 (lEttori et all I2n02h. and 



Perseus ( Churazov et al.ll2003[ ). We include all the 
observational points with two exceptions. We ne- 
glect the abundance hole observed within the cen- 
tral several kpc in Perseus (3 points at r < 10 kpc), 
Abell 262 (3 points are r < 10 kpc) and Virgo 
(1 point r < 1 kpc) . In the case of Perseus, 



Churazov et al. I (Sool) mentions that the profile 



may be affected by the bright compact source in 
this region. We also removed the data points that 
suddenly drop to very low values at large r (the 
last point in the case of Sersic 159-03 and two last 
points in the case of Abell 262). We also note 
that we used the deprojected profiles for clusters 
Perseus, Virgo, Abell 496 and Abell 4059 and pro- 
jected profiles for clusters Hydra A, Sersic 159-03, 
Abell 262 and Abell 1795 because deprojected pro- 
files were not available. We do not expect this to 
have a large effect on our conclusions, because for 
the cases for which we have both projected and 
deprojected profiles, the two types of profile are 
similar, with the main difference being that the 
deprojected profile is more noisy. The resulting 
abundance profiles are shown in Fig (21 

This figure shows that the ranges of values for 
the abundance profiles are relatively similar. We 
note that the normalization of the abundance pro- 
file of Sersic 159-03 is lower than the others. This 
lower normalization may be related to the soft X- 
ray ex cess, which is mos t likely of non-thermal 
origin (^ Werner et al. 2007) and may contribute to 
the continuum. For r < 100 kpc a central abun- 
dance peak is observed for all eight clusters while 
at larger radn (100 — 300 kpc), the profiles fiatten 
and reach an asymptotic value of a '--^ 0.2 — 0.4. 
This similar shape for all the clusters may indicate 
a common mechanism for creating the abundance 
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Fig. 2. — Observed abundance profiles in solar 
unit for our 8 clusters (Chandra or XMM). We 
decompose them into a central iron excess from the 
BCG (r < 100 kpc) and an asymptotic abundance 
from the background galaxies (r > 100 kpc). 



profiles and allows us to decompose the profiles in 
two parts: a flat contribution and a central excess. 
We note that three out of the four highest peaks 
(Virgo, Abell 4059 and AbeU 262) correspond to 
the clusters with less diffusion (see Figll] in our 
model) . 



Following the wor k of [ Bohringer et al . 
Rebusco et ail (|2005l l2006|), we subtract 



(120041 ): 

a con- 
stant value aiy from the observed abundance pro- 
files. Of, corresponds to the contribution of iron 
from the background galaxies as well from SNII. 
Our goal is to isolate the central iron excess or 
more precisely the iron contribution from the BCG 
itself. Indeed, iDe Grandi fc Molendil (|200l[ ) mea- 
sured the abundance profiles at large radii for 17 
clusters. They showed that the abundance pro- 
files of all non-cooling-core clusters are consistent 
with being constant with radius at a value of order 
0.2-0.4. On the other hand, cooling-core clusters 
present a peak near the center and an asymptotic 
value at large radii at the same value of order 0.2- 
0.4. Moreover, an an alysis of the elemental a bun- 
dance patern i n M87 (IFinoguenov et al.ll2002l) and 
Sersic 159-03 (|de Plaa et al.ll2006l) hasshown that 
SNIa dominate the metal enrichment in the central 
region and that their contribution decreases with 



radius. These observations support the idea of a 
background value from SNII in galaxies through- 
out the cluster. 

To estimate these background values af,, we 
take the average values of the last three points 
in the vicinity of 100 kpc. Then, to avoid nega- 
tive abundances, we remove the points were a < 
ttb- Since the observations of Virgo do not reach 
100 kpc and seems to be below the other curves 
at intermediate radii, we choose a typical value 
Of, = 0.2. The last observational points are lo- 
cated at a distance from the center. The values 
of ab and ri, are shown in tabled 

From the density profile and the abundance 
profile, we deduce the total observed mass of iron 
in the central excess Mpe (see table [Ij. This mass 
is the mass of iron inside rb- In our model, we al- 
ways choose svef f so that the mass of iron inside 
rb is also equal to Mpe- We also characterize the 
width of the abundance profile by computing ri/2, 
the radius which contains half the total observed 
mass of iron in the central excess (inside rb). 

5. Impact of present AGN-driven convec- 
tion on iron distribution 

Using the model described in sections [2] and [3] 
and the parameters from section [H we compute 
numerically the iron abundance profiles for the 
eight clusters in our sample and compare to ob- 
servations. For the sake of comparison, we first 
compute the abundance profile without any diffu- 
sion or inflow velocity. The result is shown as a 
dashed line in FigI3] and FigH] Clearly, such pro- 
files are too steep near the center and do not match 
the observations (diamonds). These profiles basi- 
cally show that the observed distributions of iron 
do not follow the distributions of stars within the 
BCG galaxies. This is because the effective radii of 
the luminosity distribution r^ (table [T|) are smaller 
than the observed half iron mass radii ri/2 (table 

Taking into account the diffusion and the in- 
flow velocity brings the model profiles substan- 
tially closer to the observations. This is illustrated 
by the continuous lines in FigO and Figdl In- 
terestingly, we find the same classification as in 
subsection 14.21 when studying the intensity of the 
diffusion coefficient. In the very "weak feedback 
clusters" with the least diffusion (Virgo and Abell 
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496), abundances are largely overestimated in the 
inner 20 kpc by a factor larger than four and un- 
derestimated at r > 50 kpc. "Very weak feedback 
clusters" do not have enough present convection to 
diffuse the iron to larger radii. In the "weak feed- 
back clusters" with little convection (Abell 4059 
and Abell 262), abundances are underestimated 
near 30 kpc by a reasonable factor of 1.5-2. The 
"moderate feedback clusters" with intermediate 
diffusion coefficients (Perseus and Abell 1795) pro- 
vide a much better agreement with the observed 
abundance profiles. Their present convection is 
almost enough to shape the abundance profiles so 
that they are compatible with the observational 
points. Finally, in the "strong feedback clusters" 
(Sersic 159-03 and Hydra A), the diffusion coeffi- 
cients are sufficiently large or even too large and 
the abundance near the center is compatible (Ser- 
sic 159-03) or even too small by a factor of ~ 2 
in the inner 70 kpc (Hydra A). We draw the con- 
clusion that turbulence at the level predicted by 
the AGN-driven convection model cannot be ne- 
glected when considering the iron distribution in 
"moderate" or "strong feedback clusters" . 

If AGN-driven convection is the dominant 
mechanism for transporting metals, our iron abun- 
dance profile plots in Fig [3] and Fig [4] suggest that 
Virgo- like clusters with low present convection re- 
quire a lot more convection in the past to match 
observations whereas Hydra A-like clusters with a 
lot of present convection require less past convec- 
tion to match observations. This coincides with 
the natural idea that clusters and more specifi- 
cally their central AGN were not in thei r current 
state during all their history (IPoEi|2003), an idea 
that we investigate more thoroughly in the next 
section. 



6. 



Insight into the past variations of the 
AGN povirer from the shape of the 
abundance profiles 



Pope et al. I (|2006h suggested that clusters are 



heated by Hydra-A type events interspersed be- 
tween epochs of lower AGN power. In this section 
we explore this idea and the more general possi- 
bility that the AGN power varies in some fashion 
during the lifetime of a cluster. As a consequence 
of this variation, the amplitude of the AGN-driven 
convection will also evolve. Ideally, to investigate 



the effects of a time-varying diffusion coefficient 
on the present-day abundance profile, we would 
solve equation ([6]) using the correct individual his- 
tory of convection for each cluster in our sample. 
However, we have no reliable way of determining 
these individual histories. Instead, we adopt the 
following approach. First, we develop a model for 
relating the turbulent velocity profile in a cluster 
directly to its AGN power (Sect l6.l| ). Second, for 
each cluster in our sample, we determine the con- 
stant AGN power that would be needed in order 
for the resulting turbulent velocity profile to bring 
the iron abundance profile into agreement with 
observations. We then take this constant AGN 
power to be a proxy for the time-averaged AGN 
power within the cluster over the last ~ 10 Gyr. 
In following this procedure, we make the strong 
hypothesis that the AGN-driven convection is the 
only physical process that transports iron in the 
ICM. 

One way of thinking about the cause of the 
varying accretion rate of the central AGN is pro- 
vi ded by the same AGN -driven convection model 
of lChandran fc Raseral ^07) that we have been 
using to model the turbulent velocity profile. In 
this model, the value of the AGN power (which 
we take to be comparable to Ld) is determined 
primarily by two factors: the density pcmc and 
temperature Tcoro of the intracluster plasma at 
the radius Tcorc — 50 — 100 kpc. These param- 
eters provide the outer boundary conditions for 
the strongly cooling cluster core. For larger Pcorc, 
the rate of radiative cooling within rcorc is larger, 
and thus more heating is needed to balance cool- 
ing, leading to a larger value of Lcr- (In a time- 
dependent model, if cooling were to initially ex- 
ceed the total heating, then the plasma would cool 
and fiow inward, causing the AGN mass accretion 
rate to grow until a balance between heating and 
cooling is reached.) Similarly, for a smaller Tcore, 
the conductive heating becomes smaller (the con- 
ductivity being oc T^/'^) and more convective heat- 
ing is needed in order for total heating to balance 
cooling. In the AGN-driven convection model, 
small variations in Pcove and Tcore causc large vari- 
ations in the AGN power. The value of L„ will 
thus evolve as Pcorc and Tcore change in time due 
to, e.g., mergers, streams of infalling baryons, and 
stripping of gas from cluster galaxies. 
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6.1. Analytical fit of the convection profile 
depending on the AGN power only 




Fig. 5. — Diffusion coefficient (deduced from tur- 
bulent velocity field) normalized to its maximal 
value as a function of the radius normalized to the 
radius where the diffusion coefficient is maximum. 
The dashed lines are for our 8 clusters. Their pro- 
files are self-similar. The thick continuous line is 
the fit given by Eq[8l 

In FigOwe plot -D/Dmax as a function of r/rmax 
for the eight clusters in our sample, where rmax is 
the radius at which D achieves its maximum. The 
shapes of the different curves are similar, and so we 
are able to obtain a reasonable fit to the diffusion 
coefficients using the anlalytic function 



Djr) 

Drn.n.T. 



1.6 



0.934r„ 



H -0.525 



0.934r„ 



(8) 



The result of this fitting procedure is shown as a 
solid thick line in FigO For r ^ Vrmax the right- 
hand side of equation ([8]) scales like r^^-^ whereas 
for r <C r,nax it scales like r^-^. Equation ([8]) 
matches the diffusion coefficient at r > 0.2r„iax 
fairly well. For smaller radii, Eq|5]is still compat- 
ible but the diffusion coefficient from the model 
drops more strongly for some clusters. 
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Fig. 6. — Correlation between the radius where 
the diffusion is maximum and the cosmic-ray lumi- 
nosity. Right: Correlation between the maximum 
of the diffusion coefficient and the cosmic-ray lu- 
minosity. 
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Fig. 7. — Correlation between the maximum of the 
diffusion coefficient and the cosmic-ray luminosity. 

To develop a simpler one-parameter fit for D{r), 
we note that the quantities I?max and rmax are 
closely correlated to the cosmic-ray luminosity Lcr 
of the central AGN in the AGN-driven convec- 
tion model. These correlations are shown in Fig|Bl 
FiglZl and allow us to express /?max and Tmax in 
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terms of La- using the best-fit power-law relations: 

0.44 



'^m.n.x — 3 . 1 X 



i?™„, = 2.6.10" 



10^1 erg.s-i 



kpc. (9) 



1041 erg.s-i 



0.58 



^cfo) 



The correlation between Umax and Lcr results 
from the fact that the convection is triggered by 
the gradient of cosmic-ray pressure; thus, more 
cosmic-ray energy injection will increase the con- 
vection level and the resulting diffusion. We note 
that Lcr is not a free parameter of the AGN-driven 
convection model, but rather is determined self- 
consistently within the model so that the resulting 
convective heating (which increases with increas- 
ing Lcr) is sufficiently large that the total heat- 
ing (convective plus conductive) balances radiative 
cooHng. 



10-=* 



of parameters. However, we have now achieved a 
reasonable fit of the diffusion profile (due to AGN- 
driven convection) depending only on Lcr (the 
cosmic-ray luminosity or AGN power). The evo- 
lution of the profile as a function of Lcr is shown 
FigIS] where the cosmic ray luminosity is increased 
from Virgo's value {Lcr = 5.8 x lO^^erg.s"^) to 
the Hydra A's value {Lcr = 2.4 x lO'^'^erg.s^i) by 
a factor 3.33 each time. 

The diffusion term in Eq[5] can therefore be ex- 
pressed as a function of the cosmic-ray luminosity. 
We note that the inflow term in Eql6] can also be 
evaluated from the cosmic-ray luminosity since the 
accretion rate is given by 



M = 



Lc 



0.005 X 



(11) 



We ha ve taken here the efficiency from lChandran fc Rasera 
(|2007h . This accretion rate will be useful to solve 
Eql6]in the next part. 
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Fig. 8. — Diffusion coefhcient profile (due to 
AGN-driven convection) as a function of the ra- 
dius for different cosmic-ray luminosities (AGN 
power). The cosmic-ray luminosities increases 
from 5.8 x 10"*^ erg.s^^ (Virgo on the left) to 
2.4 X lO'*^ erg.s"^ (Hydra A on the right) by a 
factor of 3.33 each time. This comes from the an- 
alytical fit EqEl and is used to determined which 
time average AGN power is required to match the 
observed abundance profiles. 

The accuracy of the analytical fit is of course 
degraded when we use equations ([9]) and (fTO]) in 
equation ([8|), since we have decreased the number 



6.2. AGN power required to match the 
abundance profile 
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Fig. 9. — Perseus abundance profile as a func- 
tion of the radius for different cosmic-ray lumi- 
nosities (AGN power). The cosmic-ray luminosi- 
ties increases from 5.8 x 10^^ erg.s"^ (Virgo on 
the middle top) to 2.4 x 10"^* erg.s"^ (Hydra A on 
the bottom right) by a factor of 3.33 each time. 
It shows that by varying Lcr we can reasonably 
match the abundance profile. 



10 



We now solve our model using the convection 
profile as given by Eql8] (for the diffusion term) and 
the accretion rate as given by Eq[Tl] (for the inflow 
term). The only free parameter is the AGN power. 
We show in FiglHlthe evolution of the abundance 
profile as a function of the AGN power. We deduce 
the luminosity L^^* which provides the best fit to 
the abundance profile. The resulting profiles are 
presented in Figl3] and Fig|4] as a dashed-dotted 
line. The first conclusion is that we obtain a good 
match for all the profiles by varying this single free 
parameter. This was not obvious because we can- 
not adjust the size and the normalization of the 
diffusion coefficient independently. Moreover, the 
diffusion drops by a factor 10 outside of the region 
0.2r„iax < r < 5r„iax and has very little effect out- 
side this region. The position of the peak of dif- 
fusion exactly where we need diffusion to explain 
the abundance profile is a striking coincidence. It 
means that by choosing the right AGN power one 
could reproduce the shape of the observed iron 
excesses. We also compute the half iron mass ra- 
dius r(y2 for each cluster (table ^ . This is the 
radius which contains half the total iron mass of 
the model inside Vb (distance to the center of the 
last observational point). The reasonable agree- 
ment with the observed half iron mass radius ri /2 
(table [2|) computed inside the same radius con- 
firms the above conclusion. 

The following question naturally arises: are 
these average cosmic-ray luminosities (Llf) real- 
istic? Nothing guarantees realistic values in our 
fitting process. The question is particularly inter- 
esting for the "very-weak feedback clusters" such 
as Virgo and Abell 262, and the "strong feedback 
cluster" Hydra A, which do not fit the observed 
abundance profiles with the present convection. 

The best-fit values of the cosmic-ray luminosity 
L^^* are summarized in table [3] and can be com- 
pared with the present cosmic-ray luminosity Lcr 
in table [1] Also shown in table [3] is the maximum 
of the diffusion coefficient D(^l^, and its radius 
'''max J which can be compared with the present 
value Tmax and Dmax in table [H We recover here 
the same kind of value for the diffusion coefficient 
value Df^l^ = lO^^ - 10^9 cm^.s"! (equivalent 
to maximum rms turbulent velocities betw een 7 

200l; 



and 110 km.s 



as 



Rebusco et all ( 2005 



Graham et al. ( 20061 ). The advantage in our ap 



co mes from the physics-based AGN-driven model 
of lChandran fc Raseral (|2007^ and that we are also 
able to derive the corresponding cosmic-ray lumi- 
nosities. 
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Fig. 10. — AGN power (cosmic-ray luminosity 
L{**) required to match the abundance profile as 
a fimction of the present AGN power (cosmic-ray 
luminosity Lcr) for our 8 clusters. The dotted line 
shows the equality between these two quantities. 
On this line, convection at present level sustained 
during the past 10 Gyr would be enough to explain 
the distribution of metals. The continuous line is 
the average < >= 4.9 x lO'*'^ erg.s""'^ over all 
the clusters. The dashe d line is the estimate from 
Voit fc Donahue (|2005 ) using entropy-floor obser- 



proach is that the shape of the diffusion coefficient 



v at ions and the dot-dashed line is an estimate from 
iBest et al. I (Hqoi) using radio observations. The 
proximity of these values to our average suggests 
an alternation of Virgo-type feedback, Perseus- 
type feedback, and Hydra-A type feedback and 
could explain the distribution of metals in clus- 
ters of galaxies. 

FigllOl shows the average cosmic-ray luminosi- 
ties L^** as a function of the present cosmic- 
ray luminosity Lcr- The dotted line shows the 
equality between these two quantities. Clusters 
located near this diagonal line have about the 
right amount of present convection to explain the 
shape of the iron distribution (if this convection 
was sustained during the past 10 Gyr). We re- 
cover the previously mentioned result that Hydra 
A has too much convection whereas Virgo and 
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A496 have too little (and the other clusters have 
approximatively the right amount of present con- 
vection). What is interesting is that the average 
< L^** >= 5 X 10'*'^ erg.s"^ is about the same 
as the average of the present cosmic-ray luminosi- 
ties < Lcr 7 X 10'*'^ erg.s"^. Moreover, the 
dispersion of cosmic-ray luminosities is much less 
since = 10^'^ — lO'''* erg.s"^ whereas Lcr = 
10'^"'^ — 10''^ erg.s^^. This is exactly what we would 
expect if in the cluster history the cosmic-ray lu- 
minosities evolved randomly between Virgo val- 
ues, Perseus values and Hydra-A values. The av- 
erage over all the clusters would be conserved and 
the dispersion of cosmic-ray luminosities would de- 
crease. 

We have searched the literature to find other 
constraints on the average cosmic-ray luminosities 
in order to check if such values are realistic. Direct 
measurement of the convection profile is, of course, 
out of reach. The better constraints come from 
feedba ck considerations. A study by IPope et al 



(l2006l) of the heating rate of cooling-core clusters 
suggests an alternation between Hydra-A type 
events and smaller-scale outflows. By studying 
the entropy floor in cooli ng-core clusters w i th an d 
without radio emission, IVoit fc Donahu^ (|2005l ): 
Donahue et al. suggest also an episodic 

heating. The invoked values are 10*^ erg.s^^ for 
10^ yr once every 10^ yr. This leads to an average 
of 10'*'' erg.s"^, which is comparable to our value 
given the uncertainties. 

The average of the observational estimates of 
the present mechanical luminosities in a large sam- 
ple of clusters could also provide a rough estimate 
of the temporal average of the AGN power (al- 
though the accuracy of such an estimate is re- 
duced by its fa ilure to take into ac c ount evolution 
with r edshift). iBirzan et aD (|2004h : lRaffertv et af 



( 20061 ) obtain a sample of 33 coohng-core clusters 
with measured mechanical luminosities. The val- 

±<J Clg.D ^ 

44 „-l 



ues range from lO"'^ erg.s ^ to 10**^ erg.s ^ with 



an average of 9 x 10 erg.s" but we note that 
the sample is certainly biased towards the high 
rate for z > 0.1. If we remove these clusters, the 
average falls to 2 x lO'** erg.s"^. These samples 
contain only a small number of clusters because 
of the difhculty of detecting X-ray cavities and 
measuring their properties. 



An other approach was taken by iBest et al 
([20071) who measured the radio luminosities for 
the BCG of 625 nearby groups and clusters. 
Using the conversion between mec hanical lumi- 
nosity and radio-lurn i nosity from iBirzan et al. 
(|2004l) : iRaffertv et all (|2006h . they converted the 



radio luminosity function into a mechanical lu- 
minosity function. The average over the sam- 
ple of BCG provides an estimate of the time- 
averaged mechanical luminosities, < Lmech >= 
2.3 x 10''2/(AL*/10" Mq) erg.s-i, with / their 
uncertainty factor, and M* the stellar mass. As- 
suming a stellar mas s to blue light ratio of 5.3 
(|Borriello et al.ll2003h . the average of our stellar 
masses is 7.4 x 10*-^ Mq. A factor f — 1 would 
mean that an energy of 1 PV is considered per 
b ubble. I f we t ake / = 7/(7 — 1) =4, as in 



Best et al.l (|2007l l. then we obtain an average me- 
chanical power of < Lmech >— 6.8 X 10'*'^ erg.s"^ 
which is in good agreement with our estimate us- 
ing the iron distribution. 



Another constraint comes from the mass of the 
central black hole, which equals the time integral 
of the mass accretion rate (plus a tiny contribu- 
tion from the initial black hole seed plus a possible 
contribution from black-hole mergers). Since the 
central black hole may have grown substantially 
before the formation of the cluster, the present 
value of the black hole mass provides only an 
upper bound on the average mass accretion rate 
of the black hole during the lifetime of the clus- 
ter. Assuming a black hole of mass 3 x 10^ Mq 
(as in Virgo), the average accretion rate must be 
smaller than 0.3 M0.yr~^. This accretion rate 
corresponds to an average cosmic-ray luminosity 
of 1.9 X 10*5(77/0.1) erg.s"*, where 77 is the effi- 
ciency with which the central AGN converts the 
rest mass energy of accreted plasma into cosmic- 
ray luminosity (which we take to be compa rable 
to the total AGN power). lAllen et all (|2006h sug- 
gest that roughly 2.2 % of the Bondi accretion 
power is transformed into AGN mechanical lumi- 
nosity. Using this value for 77, we find that the 
cosmic-ray luminosity must be less than roughly 
4.1 X 10** erg.s"*, which is larger than our value. 
Moreover, it should be noted that the Bondi accre- 
tion rate may significantly exceed the rate at which 
mass accretes onto the central black hole, since 
much of the accreting plasma may end up form- 
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ing stars (ITan fc Blackman 20051) or b eing ejected 
in an outflow ( Blandford et al.lll999f l. We thus 



conclude that our average value of is com- 
patible with present-day bla ck hole masse s . Cos - 
mological simula t ions from ISiiacki et alj pOOTI ): 



Di Matteo et al 



(2003) also indicate black hole 
accretion rates and mechanical luminosities in the 
same range of values with large variations. In 
these simulations, the final black hole mass comes 
from mergers and accretion in high redshift quasar 
phases. 

Our value of L^^* ^ 5 x lO'*'^ erg.s^^ seems rea- 
sonable and compatible with current constraints 
(see Fig flU)) . If our scenario is correct, it indi- 
cates that the iron distribution could put strong 
constraints on the AGN feedback model. For 
instance, an average AGN power greater than 
2.4 X 10**^ erg.s"^ would destroy the observed 
abundance profile of most of the clusters as illus- 
trated FiglQl Moreover, our scenario gives fur- 
ther support for the idea that convection plays 
a fundamental role in clusters of galaxies by ex- 
plaining the density and temperature profile as 
well as the shape of the iron abundance excess. 
Our model is, however, very simplistic and re- 
lies on many approximations. In order to be 
more realistic, one would have to carry out cos- 
mological simulations with AGN feedb ack and 
supernovae enrichment as Sijacki et al.l ([2007); 
Di Matteo et al.l (|2007l ). However, a proper treat- 
ment of the convection would also require the 
inclus ion of cosmic rays and a nisotropic trans- 
port ( Rasera fc Chandran 2007h since these two 
ingredients modify the convec tive insta bility crl 
terion fea lbus 2000; Chand ran fc Denn is 20(>_ 
Dennis fc ChandranI boOSi : Parrish fc Ouataer' 
20081: lQuataertll2008l ) 



7. Conclusion 

In this article, we have studied the impact 
of AGN-driven convection on the shape of the 
iron abundance profile for 8 cooling-core clus- 
ters: Perseus, Hydra A, Sersic 159-03, Abell 
262, Abell 1795, Virgo, Abell 496 and AbeU 
4059. We have used th e iron injection model from 
Bohringer et al.l (|200J) where metal production by 
SNla and winds follows the Hernquist light pro- 
file of the central brightest cluster galaxy (BCG). 



We have al so used the steady- s tate convection 
model from Chandran fc Raseral (|2007t ) to deter- 
mine the turbulent velocity profile in the ICM. In 
this model, AGN-driven convection is the domi- 
nant mechanism for transferring energy from the 
central AGN to the ICM thereby preventing large 
quantities of plasma from cooling to low temper- 
atures. 

Stochastic motions with an rms radial velocity 
fluctuation unl correspond to a diffusion coeffi- 
cient D = 0.5l uNr, in the mixing lengt h-theory 



employed bv IChandran fc Raseral ( 2007f ). where 
I is the mixing-length which we set to 0.4r. We 
have solved a ID advection-diffusion equation for 
the iron abundance profile, and compared our re- 
sults to XMM and Chandra abundance profiles. 
The profiles obtained without diffusion (metal in- 
jection only) are too peaked toward the center. 
Taking into account AGN-driven convection im- 
proves the abundance profile because it smoothes 
the center. For most objects in our sample the 
modeled profiles do not differ greatly from the ob- 
servations. However, the less convective clusters 
(Abell 496 and Virgo) require much more convec- 
tion in order to match observations whereas the 
most convective cluster (Hydra A) requires less. 

Making use of the approximate self-similarity of 
the diffusion coefficient profiles in the AGN-driven 
convection model, we model the diffusion coeffi- 
cient as a function of radius r that depends only on 
one parameter, namely the cosmic-ray luminosity 
(AGN power). We have used this function to find 
if any reasonable AGN power could provide a good 
fit to the observed iron distributions. Although 
the diffusion coefficient is non-negligible only in a 
limited range of radii around its maximum, and 
although we cannot adjust the amplitude and the 
position of this maximum independently, we find 
a good match for all the clusters. 

The best-fit AGN powers i^^* can be thought as 
the time average over the cluster history (10 Gyr) 
required to explain entirely the width of the abun- 
dance profile. We found a range of values for 
between 10**^ and 2x 10^"* erg.s"^, with an average 
of 5 X 10"*^ erg.s"^. Such a value seems quite rea- 
sonable and compatible with current constraints 
on the entropy profile, black hole mass, and radio- 
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inferred average mechanical luminosity. The shape 
of the abundance profile could therefore be a fossil 
record of the past AGN-driven convection. 

This model relies on several approximations. 
The separation of the abundance profile into a con- 
stant component ai, (from the background galax- 
ies and SNII) and an iron excess (from the central 
BCG) is one of the major approximations. A bet- 
ter method would be to also include the contribu- 
tion from other galaxies and SNII and to compare 
with the total abundance profile, however this con- 
tribution also suffers from a lot of uncertainties. 
We also neglect all other possible sources of con- 
vection such as supernova winds, the jet itself, and 
stirring by infalling galaxies. This may explain 
the discrepancy between our results and the ob- 
served abundances at the very centers of several 
of the clusters, in some of which an abundance 
hole is observed. Another step towards a more 
realistic model would be cluster simulations with 



meta l injection and AGN feedback (jSiiacki et al 



20071 : iRoediger et al.ll2007[ ) but also cosmic-ray in 



jection, anisotropic transport and cosmic-ray dif- 
fusion. By modifying the convcctive instability 
criterion these three last ingredients should drive 
turbulence to a level close to the one predicted by 
the mixing length theory. 
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Fig. 3. — Abundance profile in solar unit as a function of the radius in kpc for Virgo, Abell 262, Sersic 159-03, 
and AbcU 4059. The dashed hne is the resulting abundance profile without diffusion. The continuous line 
is the resulting profile with the present cluster diffusion. The dot-dashed line is the profile with the best fit 
cosmic-ray luminosity value and subsequent diffusion. The diamonds are the observations. Past and present 
AGN-driven convection seems to play an important role for the transport of iron. 
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Fig. 4. — Abundance profile in solar unit as a function of the radius in kpc for Hydra A, Abell 496, Abell 
1795, and Perseus. The dashed line is the resulting abundance profile without diffusion. The continuous line 
is the resulting profile with the present cluster diffusion. The dot-dashed line is the profile with the best fit 
cosmic-ray luminosity value and subsequent diffusion. The diamonds are the observations. Past and present 
AGN-driven convection seems to play an important role for the transport of iron. 
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Cluster 


z 


(h-iMo) 


(erg.s"^) 


BCG 


Lb 
(L|) 


r 

(kpc) 


Virgo 


0.00372 


3.1 X lO^'' 


5.8 X 10"! 


M87 


6.4 X lO^*^ 


6.1 


Abell262 


0.0155 


3.2 X 10^^ 


6.6 X 10"^ 


NGC708 


3.8 X 10^" 


24 


Sersicl59-03 


0.0572 


9.6 X 10^^ 


1.1 X lO'*'^ 


ESO291-009 


1.3 X 10" 


27 


Abell4059 


0.0466 


6.8 X 10^"^ 


1.7 X 10"-' 


ESO349-010 


1.9 X 10" 


22 


HydraA 


0.0550 


7.1 X 10^^ 


2.4 X lO'''^ 


3C218 


2.6 X 10" 


44 


Abell496 


0.0322 


8.2 X 10^^ 


1.9 X 10"*^ 


PGC015524 


1.6 X 10" 


50 


Abelll795 


0.0639 


5.0 X 10^* 


1.0 X 10''* 


PGC049005 


1.1 X 10" 


42 


Perseus 


0.0179 


6.2 X 10^"^ 


9.8 X 10*-' 


NGC1275 


1.7 X 10" 


15 



Table 1: Clusters and BCGs general properties 



Note. — z is the redshift, M^ir is the virial mass, Lcr is the cosmic-ray luminosity (AGN power), BCG is the brightest 
cluster galaxy. Lb is the blue luminosity of the BCG and re is the effective radius which contains half the projected luminosity. 
Methods, definitions and references are described in section |4] 



Cluster 


^max 

(kpc) 


I-^ max 

(cm^.s^-') 


n 
(kpc) 


ab 

(a©) 


Mpe 
(Mq) 


ri/2 
(kpc) 


Virgo 


4.4 


8.0 X 10^^ 


58 


0.20 


4.2 X lO-" 


36 


Abell262 


62 


6.3 X 10^« 


200 


0.29 


1.4 X 10** 


55 


Sersicl 59-03 


98 


2.6 X lO^'-* 


260 


0.14 


2.4 X 10» 


110 


Abcll4059 


39 


4.7 X 10^** 


85 


0.43 


1.9 X 10** 


44 


HydraA 


78 


2.5 X 10^'-* 


140 


0.28 


2.0 X 10** 


57 


Abell496 


7.8 


8.0 X 10^"^ 


200 


0.28 


2.6 X 10** 


80 


Abelll795 


67 


9.9 X 10^** 


170 


0.33 


4.8 X 10** 


92 


Perseus 


30 


9.4 X 10^** 


98 


0.46 


1.7 X 10** 


53 



Table 2: Diffusion properties derived in the AGN-convection model and abundance profile parameters 



Note. — rmax is the position of the maximum of the diffusion coefficient and Dmax is the maximum, ri, is the radius of the last 
bin of the observed abundance profile, ai, is the background abundance (relative to the solar abundance of lAnders &: Grevessj 
dlQSgD ) which is subtracted to the abundance profile to obtain the iron excess, Mpg is the total observed iron mass, and ri/2 is 
the iron half mass radius (which therefore contains an iron mass of M^e/2). Methods, definitions and references are described 
in section |4] 
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Table 3: AGN power and diffusion properties for the best fit model to the abundance profile 



Note.— lS* is the cosmic-ray luminosity (AGN power) required to reproduce the observed abundance profile. From this 
AGN power, we derived the diffusion coefficient profile, r^"^^ is the radius of the maximum of the diffusion coefficient profile 

and f is the maximum, sr^ff is the effective SNIa rate (which includes the wind contribution to the iron injection), r^^^* is 
the iron half mass radius of the model (which contains half the total iron mass inside r^). Methods and definitions are described 
in section |2] and section |6] 
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